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Abstract. We show that an approximate solution to the amended non-linear Balitsky-Kovchegov evolution 
equation which was formulated for hard (large Q 2 ) QCD processes, can be extended to provide a good 
description of photoproduction and soft hadronic (non perturbative) reactions. 
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1 Introduction 

Further proof that the linear QCD evolution equations fail 
to describe the DIS data for low values of Q 2 has recently 
been provided by the ZEUS collaboration P^. See Fig.l. 
For values of Q 2 > 1 GeV 2 the ZEUS NLO QCD fit pro- 
vides an excellent description of the data, however, as Q 2 
becomes smaller the discrepancy between the predictions 
and the data increases. 

The reason for this is well known, and is due to the fact 
that in the linear evolution equations DGLAP and BFKL, 
the splitting functions only incorporate the production of 
partons (gluons) . GLR [2] in their classical paper pointed 
out that when the density of partons becomes sufficiently 
large (that they overlap), one has to include non-linear 
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(annihilation) processes in the evolution equations. A re- 
cent application of the ideas of GLR using the dipole for- 
malism has been suggested by Balitsky [3] and Kovchegov 
0] (which we will denote by BK). The advantage of the 
BK equation is that: 

(i) it accounts for saturation effects due to high parton 
densities; 

(ii) it sums higher twist contributions; 

(iii) it allows one to extrapolate to small values of Q 2 
(large distances). 



2 Amended Non-linear BK Equation 

The BK equation can be written in the form ,5; 
dA>(x iA Y) 



C F a s 



dY 



ip ^-02^-12 

• ( 2A>(x 02 ,6,r) - N( X02 ,b,Y) N( Xl2 ,b,Y) ) (1) 

where N(r±, b, Y) denotes the imaginary part of the am- 
plitude of a dipole of size r± elastically scattered at impact 
parameter b, i.e. 



N(r ± ,b,Y) = Imai pole (r ± ,b,Y). 



(2) 



Fig. 1. F 2 data at very low Q 2 compared to ZEUS-S NLO 
QCD fit, from ref.[l]. 



Y = Iti(1/xbj) denotes the rapidity. In Eq. (1) the linear 
term corresponds to the LO BFKL evolution, while the 
non-linear negative term is responsible for unitarization. 

A deficiency of the BK equation is that it does not con- 
tain evolution in Q 2 , and therefore lacks the correct short 
distance behaviour. To remedy this we have introduced a 
correcting function AN(r±, b, Y) which accounts for the 
DGLAP evolution in Q 2 . Our full solution therefore con- 
sists of the sum of two terms 

N(r±,b,Y) - N(v ± ,b,Y) + AN(r ± ,b,Y) (3) 
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Full details of the calculation are contained in ref. 
where we have simplified the treatment of the dependence 
on the impact parameter b. We solve for the initial condi- 
tion at b = 0, and then restore b dependence by assuming 
factorization and that the dipole profile function inside the 

_ i> 2 

target is given by a Gaussian distribution S(b) = e ~r? for 
all b. This is obviously an over simplification (see discus- 
sion in Section 5), and the true impact parameter depen- 
dence of N(r±, b, Y) is far more complex. 

The deep inelastic structure function F 2 is related to 
the dipole cross section (defined in Eq.(2)) 



F 2 (x BJ ,Q 2 



d 2 r t 
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At: 2 

CTdipole(r±,XBj) 



dzP^(Q 2 ;r ± ,z) 



(4) 



Where P 1 (Q 2 ;r±,z) denotes the probability of the de- 
cay of a virtual photon having four momentum Q 2 into 
a colourless dipole (qq pair) of size r_i, with the quark 
(anti-quark) taking a fraction z (1-z) of the virtual pho- 
ton's momentum. 



3 Numerical solution of the equation for DIS 

In making our fit we include all available data satisfying 
the following criteria 10~ 7 < xbj < xqbj = 10 -2 and 
Q 2 > 0.04 GeV 2 . We solve Eq.(l) as an evolution equa- 
tion in rapidity with a fixed grid in rj_ (= 2/Q GeV) space 
and a dynamical step in rapidity. We fit to 345 data points 
for F 2 (xbj, Q 2 ), and obtain an overall x 2 /^/ ~ !• I n re f- 
[HJ we also compare to data on the logarithmic slopes of 
F2 he. j} F A± and ,?! n f? s . The values that we obtained 

z dh\Q z d(lnl/x) 

for A = ^/"/L were A w 0.07 for Q 2 = 0.05 GeV 2 and 

0.3 for Q 2 = 150 GeV 2 



d{lnl/x) 

xbj = 10~ 5 , increasing to A 

and xbj = 4.10" 2 . i.e. our formalism was successful in 
describing not only the short range (large Q 2 ) data but 
also the "soft" data ( at very low Q 2 ) where the tradi- 
tional Pomeron intercept is ~ 0.08. 



4 Extension of BK formalism to 
photoproduction 

The surprising results discussed above, that with the 
amended BK equation we found an excellent description 
of all DIS data for 0.05 < Q 2 < 200 GeV 2 prompted the 
question, whether this formalism could also successfully 
describe soft processes e.g. photoproduction? 

To extend our formalism to photoproduction 6 , it is 
necessary to make the following alterations: 

(i) We need to introduce a finite mass as a cutoff for the 
r± integration in Eq.(4 ), we take this parameter as m q 
(the quark mass). 

(ii) The variable x (= Xbj) is not defined for 7-p scatter- 
ing, and we relate x to the energy by introducing a non- 

perturbative scale Qq, and taking x — ■ To reduce 
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Fig. 2. 7-p and DIS cross-sections at very low Q 2 . Solid line 
VQ and dotted line SR parametrization for low energy contri- 
bution (see text for details). 



the number of free parameters we have set Qq = 4 m 2 . 
On fitting to the high energy photoproduction data we 
find a value of m q — 0.15 GeV. In the colour dipole for- 
malism one can only hope to reproduce the asymptotic en- 
ergy dependence i.e. the Pomeron contribution. We also 
need to include a non-singlet (NS) term to account for 
lower energy (higher x) contributions. We test two possi- 
ble forms for the NS term: 
1) Based on Valence Quark Model (VQ): 



F, 



VQ 



■ Q 2 s(i +0 ) (1 GeV 2 + /x 2 ) 
I GeV 2 ' (Q 2 + /i 2 ) 

■ ]T e 2 qY(Q 2 = I GeV 2 ) 



(5) 



i—u.d 



We freeze the CTEQ6 valence quark contributions at Q 2 = 
1 GeV 2 , and on fitting to data, the best fit parameters are 
(3 = and fi 2 = 0.13 GeV 2 . 

2) Based on exchange of Secondary Reggeons (SR): 



'SR 



/(0) 



M 2 



(Q 2 



M 2 



{-T R 

so 



with a R — —0.45 and s = 1 GeV 2 , /(0) denotes the 
residue at Q 2 — 0. Fit results in /(0) = 0.19 mb and 
M 2 = 2 GeV 2 . 

The results of the fit to photoproduction and low Q 2 
DIS data with the two alternate low energy contributions 
are shown in Fig. 2. 



5 The BK formalism and soft hadronic 
processes 

Based on the successful extension of our BK formalism to 
photoproduction (Sec. 4), we also applied the procedure to 
soft hadronic interactions [Jj- This necessitated adapting 
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Fig. 3. K p and 7r~p cross-sections. The full line is the predic- 
tion in our model and the dotted line using the Golec-Biernat 
Witsthoff dipole parametrization |10) 
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Fig. 4. pp and pp total cross-sections. The full lines are the 
model predictions with a Gaussian profile for the impact pa- 
rameter dependence. The dashed line for a K\ profile and dot- 
ted line for the GBW dipole model. 



the basic formula Eq. (4) of the dipole model, and hypoth- 
esizing that the hadron-proton cross-section is given by: 

Vhadron-proton(x) = J d 2 r± | 1ph(r±) | 2 <Jdipole(r X , x) 

where tph(j'±) represents the wave function of the hadron 
which scatters off the target proton, and the energy de- 

pendence is given by x = , with Q\ being an additional 
non-perturbative scale. 

For the hadronic wave functions we use the form sug- 
gested by the Heidelberg group, Dosch et al. jHJ. The 
hadronic transverse wave function is taken as a simple 
Gaussian, where the square of the wave function is given 

by 

I ^(rx) | 2 = -4r«sp(-^) (7) 

and Sm is a parameter related to the meson size. We have 
used S-n — 1.08 fm and Sk = 0.95 fm, which is related 
to the electromagnetic radii. The exotic channel K + p 
has no secondary Regge contributions, while for the tt~ 
p reaction we have added a secondary contribution a la 
Donnachie and Landshoff |5] . Our results are displayed in 
Fig. 3, for more details regarding the fit shown in Fig. 3 we 
refer the reader to 0. 

For a baryon projectile, we assume that the baryon 
is constituted of two colour dipoles, one dipole formed 
around two quarks, and the second dipole from the cen- 
tre of mass of these two quarks to the third quark in the 
baryon. Our results for p-p and p-p scattering are shown 
in Fig. 4. The assumed Gaussian dependence in b, corre- 
sponds to a e _t behaviour of the differential cross-section. 
The full lines in Fig. 4 show that this is a poor description, 
a much better fit is obtained by assuming that the momen- 
tum transfer dependence is of a "form factor" dipole type, 
which transforms to a Bessel function K\ dependence in 
impact parameter space. See ref. jj] for more details. 



6 Conclusions 

We have shown that our approximate solution to the amended 
non-linear BK equation which was formulated for hard 
(large Q 2 ) processes [5], can be successfully extended to 
describe photoproduction and soft (non perturbative) hadronic 
reactions. Two outstanding problems remain: 

(i) the impact parameter dependence of the colour dipole. 
The assumption of factorization in b space, plus impos- 
ing a Gaussian like behaviour in b (for all values of b), 
is obviously naive, as can be judged from the results we 
have obtained for pp scattering. The search for the correct 
impact parameter dependence of the solution to the BK 
equations continues. 

(ii) the form of the hadronic QCD wave function. The 
form suggested by [5J, is only a first approximation and 
should be improved. 

I would like to thank my colleagues Jochen Bartels, 
Genya Levin, Michael Lublinsky and Uri Maor for a most 
enjoyable collaboration, the fruits of which are presented 
here. 
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